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ABSTRACT 

Rotating deformed neutron stars are important potential sources for groundbased 
gravitational- wave interferometers such as LIGO, GE0600 and VIRGO. One mecha- 
nism that may lead to significant non-asymmetries is the internal magnetic field. It is 
well known that a magnetic star will not be spherical and, if the magnetic axis is not 
aligned with the spin axis, the deformation will lead to the emission of gravitational 
waves. The aim of this paper is to develop a formalism that would allow us to model 
magnetically deformed stars, using both realistic equations of state and field config- 
urations. As a first step, we consider a set of simplified model problems. Focusing 
on dipolar fields, we determine the internal magnetic field which is consistent with a 
given neutron star model. We then calculate the associated deformation. We conclude 
by discussing the relevance of our results for current gravitational- wave detectors and 
future prospects. 



1 INTRODUCTION 

It is well known that a non-axisymmetric deformation of a rotating neutron star will lead to a time varying quadrupole moment 
and could provide a good source of gravitational waves. In fact, rapidly rotating neutron stars are important potential sources 
of continuous gravitational waves for interferometric detectors such as LIGO, GEO600 and VIRGO as well as the planned 
next generation interferometers which should be able to target them specifically. There have even been su ggestions of narrow- 
bandi ng Advanced LIGO in order to target the Low Mass X-Ray Binaries (LMXBs) [see for example iBradv fc Creighton 
a s it is thought that gr avitational waves may be playing a role in setting the spin equilibrium period of these systems 



Bildstenlfliii : Undersson etaL2005. ) 



There are a number of mechanisms that may lead to a neutron star being deformed away from symmetry. First of all, the 
crust of a neutron star is elastic and can support "mountains". The size of the deformation that can be sustained depends 
on many factors, such as t he e quation of state and the evolutionary history of the crust. This problem has been studied 
bv lUshomirskv et al.l (|20od ) and lHaskell et al.l (|2006h . Another mechanism for producing asymmetries is an oscill ation mode 
develo ping an instabili t y, dri ven by gravitational radiation reaction, such as the r-mode instability proposed by AnderssonI 
1 1998t ) (see Andersson 1 19981 ) for a relatively recent review). Deformations can also be caused by the magnetic field. Neu- 
tron stars are known to have s ignificant magnetic fields and it is well known that a magnetic star will not remain spherical 
I Chandrasekhar fc Fermi 19531 ). If the magnetic axis is not aligned with t he rotation axis the deformation will not be ax- 
isymmetric and this will lead to gravitational wave emission. ICutlej 1I2OO2I ) has suggested that a strong toroidal field could 
force a precessing neutron star to become un stable and " fiip" to become an orthogonal rotator. This would be an optimal 
configuration for gravitational-wave emission. iHevll (j2000l ) has suggested that magnetised white dwarfs may be interesting 
sources of gravitational waves. Several numerical st udies have been directed at unders tanding t he gravitational-wave em ission 
of magnetically distorted neutron stars, for example iBonazzola fc Gourgoulhon j 19961 ). Finally, iMelatos fc Pavnd (|2005l ) have 
studied the closely related problem of magnetically confined accretion, which may lead to large deformations. 

The aim of this paper is to investigate deformations due to the interior magnetic field in more detail. We want to develop 
a framework that would allow us to quantify the relevance of these asymmetries for any given stellar model and magnetic 
field configuration. As a starting point it makes sense to study some simple model problems. This will give us a better idea of 
the nature of the problem and what the key issues are. Moreover , it is a natural way to proc eed given that the interior field 
structure is uncertain. The best current models [see, for example, Braithwaite fc Spruit 1 20061 )] tell us that the field will tend 
to have a mixed poloidal and toroidal nature. Any analysis should be able consider such generic configurations. 
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It is, of course, the case that magnetically deformed stellar models have been studied previously and we can benefit greatly 
from the existing literature. It is particularly important to appreciate that the magnetic field configuration is constrained for 
any given stellar model. In effect, the magnetic field must be solved for together with the fluid configuration. This statement 
is quite obvious, but it has not always been appreciated in discussions of gravitational waves from magnetically deformed 
stars. Hence, we feel that it is appropriate to pay special attention to it here. This means that our paper has two main parts. 
The flrst part, which comprises Sections 2-4, is focused on the issue of the permissible field configuration. In many respects 
it is an adaptation of already existing results. It is nevertheless important to have this discussion, since it provides the main 
input for the deformation calculation. This calculation, the second part of the paper, is presented in Section 5. The potential 
impact of our results on gravitational-wave observations is discussed in the concluding section. Throughout the paper, we 
consider only dipolar fields, but the developed formalism is general and can easily be extended to any magnetic field. 



2 MAGNETIC FIELDS IN STELLAR INTERIORS 

Chandrasekhar fc Ferniil (|l953l 'l where the first to realise that a star would not remain spherical in the presence of a strong 



magnetic field. They calculated the deformation of an incompressible star with a constant dipolar field by minimising the 
energy of the configurat ion. The case o f a c onstant density star with an internal poloidal field matched to an external dipole 
was later considered by Ferrarol (1954) and GoosensI (1972), by solving the Euler equations. In our analysis we shall use the 



latter approach. However, before considering magnetic deformations it is important to understand which are the permissible 
field configurations, given the equation of state. We shall see that the range of permitted fields is quite restricted, a fact 
which is often overlooked when discussing magnetic deformations. Let us, therefore, move on to describing the equilibrium 
configuration of a magnetic star. We will do this by assuming that the magnetic energy is small compared to the gravitational 
energy and that magnetic effects can be treated as a perturbation of a spherical, non magnetic, background. This should 
always be the case. 

The equations of hydrostatic equilibrium are: 

Vj5 „^ (V X B) X B L 

— + = = - — (1) 

where B is the magnetic field, p the pressure, p the density and L defines the Lorentz force. The gravitational potential "I> 
obeys Poisson's equation 

= iivGp. (2) 
The magnetic field must also obey, from Maxwell's equations, 

V ■ B = (3) 
Following Roxburgh 1 196d ) we take the curl of equation ([T|. This leaves us with an equation for the magnetic field 

„ /Bx(VxB)\ „ 

V X =0 (4) 



P 

Since this equation contains the density it should be solved simultaneously with equation In other words, the magnetic 
field structure is constrained by the density profile. 

We now assume that the magnetic field only produces small deviations from a spherically symmetric background model 
(essentially, we assume that the ratio of magnetic to gravitational potential energy is small) . This allows us to expand all our 
variables in the form 

iP{r,e)^Mr) + Mr)Pi (5) 

where Pi are the standard Legendre polynomials and tpi is a small perturbation of 0{B^). In the following we will concentrate 
on quadrupole (I = 2) deformations, simply because they are optimal from the gravitational-wave emission point of view. 
However, the formalism applies to the general problem, in which case the perturbation is given by a sum of Legendre 
polynomials. 

We can first of all solve the structure equations in the absence of a magnetic field and obtain the background model. The 
result is then fed into equation Q to determine B, which we will need to solve equation ([1]) and ((2]) to first perturbative order 
for the quantities p,p and Restricting ourselves to axisymmetry, the (f) component of the magnetic force must be zero, as 
there is nothing to balance it in equation ([l}. Hence 

[B X (V X B)]^ = (6) 

With this understanding, let us examine some general magnetic field solutions for various background models. If one splits 
the magnetic field into two components, a poloidal one Bp = {Br, Eg, 0) and a toroidal one Bt = (0, 0, B^), and introduces a 
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stream function S(r, 0) such that 
1 dS 



Be = 



■sine de 
1 dS 



rsinO dr 

then equations (|3]) and ([6| reduce to ( Roxburgbl 1966l l: 

Bp • V(rsin6IBt) = 

which gives 

R /3(5) 



(7) 
(8) 
(9) 



r sin e 

Where /3 is some function of the stream function 5*. This means that the toroidal part of the field is a function of the poloidal 
part. Equation Q then takes the form 



dS 



d 



dr \ pr sin 9 86 



1 d'S Id 



-sin^ dr^ 



1 d 



do \ sin 9 



1 dS 



P dp 



pr-' sm 



1 dS 

sin e de 

= 



+ 



P dp 
pr'^ sin^ Q de 



d 



dS 



pr sin 9 dr 



1 d^s 

■sin^ dr^ 



(10) 



In the following we shall focus on dipole solutions. We can then take a stream function of form 

S {r, 9) = A(r) 8111^9 (11) 

or directly look for solutions of the form: 

B = r{W{r) cos 9} + 0{X{r) sm9} + 4>{iZ{r) sm9} (12) 

where the phase of the (f> component is chosen in such a way as to produce real valued equations in the following. For this 
kind of field, equation (|3} gives 

rW'{r)+2[W{r)+X{r)] (13) 

where the prime indicates differentiation with respect to r. We now need to solve equation (|4]) [or equivalently equation (|10p ]. 
which depends on the density profile p{r). It is thus necessary to prescribe an equation of state for the stellar matter. 



3 UNIFORM DENSITY STARS 

In order to begin investigating the problem and understand how restricted our choice of magnetic field configurations may 
be, it is useful to consider a constant density star. Insert p=constant into Q to get 

W {rZ) ' + 2ZX = (14) 
W{r'^Z" -2Z)->r2Z{rX' -X)=Q (15) 
W[r'^X" - 4:{W + X)\+2Z{rZ' - Z) =Q (16) 



3.1 Poloidal fields 

Let us first restrict our attention to purely poloidal fields. In this case we can combine equation (|16p with equation ([3]) to get 

r'^X'" + 4rX" - 4X' = (17) 

which has the solution 

X = D + Cr^ (18) 

with C and D constants and where the oc r~^ solution has been excluded to ensure regularity at the centre. This gives, for 
W, 

W = ~D~ ^Cr^ (19) 

Note that for C — this corresponds to a uniform dipole field. The same solution can be obtained by solving for the stream 
function. It is then sufficient to take S{r, 9) of the form in equation (1111) and take /3 = so that the magnetic field reduces to 

„ /2Acos6' -A' sin 61 
B = 5 , , 



(20) 
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We can now solve equation (|10p with the boundary conditions that the field must remain finite at the centre, i.e. 

A A' 

— finite at r = 
r-' r 

while at the surface the field must be continuous with an external curl free dipole field, so that we have 

\-A' = , at r ^ R. 

r 



(21) 



(22) 



The solution to this problem was given by iFerrard ([195J), who considered a constant density star with a current density of 
the form 



where B is a constant parameter. In this case equation (|10|) reduces to 

V^S=^r%in^( 



which yields for the stream function: 



S = B 



3 



(23) 



(24) 



(25) 



The magnetic field is then 

Br — B 



5 \R/ 



1 1 / r 
3 ~ 5 

3 5 \RJ 



Be = B 
which corresponds to taking 



^ = -5^ D = - 



(26) 



(27) 



in the previous solution P8|) -(I19 | ). 

In Section 5 we will work out how this field deforms the shape of the star. As we will consider several other field 
configurations, and it would be useful to be able to compare the results, it is worth thinking about how such a comparison 
would be carried out. After all, each model field will be naturally described by some set of parameters which may not be easily 
translatable from model to model. Intuitively, one would e xpect the deformation of the star to depend on the ratio of magnetic 
to gravitational potential energy [see for example ICutlerl (2003)]. Hence it would make sense to consider the magnetic energy 
as a useful measure. Introducing the volume averaged magnetic field energy B, one can show that for the field in (f26)l we have 
B = O.IS. 



3.2 Toroidal field 

For a purely toroidal field, on the other hand, the equations reduce to 

rZ' = Z (28) 

which gives 

Z = ~icr (29) 
This corresponds to a uniform current distribution inside the star, i.e. |V x B| = 2c =constant, pointing along the z-axis . 



3.3 Mixed poloidal/toroidal field 

Let us return to the general case. If we formally solve equation (I13p for X, then equation (|14|l becomes 



W 



W 



Z = arW 



Hence, unless W — or Z = we have 



where a is a constant. Note that equation 1)15^ is implied by (|14|l by the use of (|13|) : thus the substitutions 

Z^arW and X = ~W ~ ^w' 
solve equations (|14|l and (|15|l . The last equation, (I16|l . then becomes 

^rW (r^w'" + 4rW" - r(l + aV^)w') = 



(30) 

(31) 
(32) 

(33) 
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Thus we find that the solution can be written 

= ^[Ci(2ar-l)e''"- + C2(2ar + l)e-'n+C3 (34) 

leading to 

X = --!--[Ci(4aV' - 2ar + l)e''"' - C2(4aV^ + 2ar + l)e-''"'] - Cs (35) 

where d are constants. In order to ensure regularity at the centre we must take C'l ~ C2, and the central values of the fields 
(which we shall denote with a subscript c) then become 

Wc = -Xc = y C-ia^ + C3 and = (36) 

As we can see, the parameter a is just a length scale and can be absorbed in the other constants if we define a new dimensionless 
radial coordinate x = 2ar. This allows us to redefine the constant 

Ci ^ Ci/(2a)' (37) 

The free parameters in the regular solution now correspond to the new Ci and C3 and are only two (if we exclude the trivial 
choice of scale or units given by a). Explicitly we have 

2Ci cosh X 2Ci sinh x 
W — h 63 

X-' X'^ 

Cicosha; Ci (a::^ + 1) sinh a; 
X = G3 

X'^ x"^ 

CiCosha; Cisinha; 1 

Z = h 7:2:03 (38) 

X x-^ 2 

We can further interpret the parameters by noting that the central values of the fields are Wc = —Xc = |Ci + C3 {Zc ~ 
still). We can thus use Ci as an overall scale for the field and define W — W/Ci and likewise for the other variables. Thus, 
using C3/C1 = Wc — 2/3, we obtain 

W = ^^^^-^ + Wc~2/3 

x'' X'' 

X = £2!^_(^!±4£^-w., + 2/3 

' cosh a; sinh a; 1 , , , , 

^ = — ^+2^(^^-2/3) (39) 

We now have only one non-trivial parameter, Wc, which controls the ratio of poloidal to toroidal field, while the other 
parameters a and Ci have been absorbed into the definition of our variables. They can simply be interpreted as scales of the 
problem (specifically a length scale and the scale of the magnetic field). 



3.4 Boundary conditions 

Let us now focus on the boundary conditions that the magnetic field needs to satisfy at the surface of the star (at the centre 
it is sufficient to impose regularity). The solenoidal nature of the field requires continuity of the radial component 

(B") = (40) 

where (B^) = BJ^t ~ A'lit indicates the discontinuity between the external and internal parts of the field at the surface. 
Furthermore, we shall require continuity across the boundary of the traction vector = T^-'fij, i.e. the projection of the total 
stress tensor along the normal unit vectofl This is equivalent to requiring local force balance per unit area. Our condition is 
hence 

it') = (41) 
The stress tensor is the sum of a fiuid piece and a magnetic piece 

T" = -pS'' + ^(b'B' - Ib^S''] (42) 
47r y 2 J 

Let us consider the projection along the normal to the surface. The normal vector to the perturbed surface will have the form 

ns =r + ns (43) 



^ While the main part of our discussion uses standard spherical coordinates, the discussion of the boundary conditions becomes clearer 
if we use a coordinate basis. This means that the Einstein convention of summation over repeated indices applies. 
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where fig indicates the correction of 0{B^). Projecting T'^ along this vector, we obtain 



Atv \ 2 J 



(44) 



Where p'^ is the background pressure and 5p the first order perturbation. Note that the magnetic term is already 0{B^) and 
one can take fis ~ r for this term. At the surface the background pressure p" vanishes, so we must impose continuity for 

= -Sp+l-[{B''f~^B'] (45) 

t" = —B'B" (46) 
47r 

t'i' = —B'B''' (47) 
47r 

As we aheady have {_B' ) — 0, the 9 and (f) components of equations (I47|l demand that (B^) — {B'^) — 0, which then leads to 
(B^) = 0. As p = in the exterior, it must be that Sp = at the surface (all quantities are now 0{B^) so we can consider the 
surface of the unperturbed configuration). The conclusion is thus that all components of the magnetic field must be continuous 
across the interface, i.e. that we have no surface currents. 

We should stress that we are in principle allowed to introduce surface currents at this point. This would lead to disconti- 
nuities in the 9 and (p components of the field. However, we feel that unless there are physical arguments dictating the nature 
of such currents there is no reason to introduce them. If we allow for the presence of an arbitrary current at the surface we 
have too much speculative freedom and it is not clear to what extent various models make sense. Hence, we find it more 
natural to restrict ourselves to models where the above traction conditions apply. 

Note that the traction conditions are automatically satisfied if we have a purely toroidal field. The external magnetic 
field, in fact, then solves 

V-B = VxB = (48) 

and the assumption of axisymmetry forces a vanishing toroidal component B*^ — 0. In the case of a purely poloidal field it is 
sufficient to match with an external curl-free dipole. In the case of a mixed poloidal/toroidal field, however, one must have 
that at the surface 

Z(R) = (49) 

and as we have found that Z — arW , cf. (|31[). it must also be the case (provided that a 7^ 0) that W = G. This means that, 
under axisymmetry, a mixed poloidal and toroidal dipolar internal field and a general multipolar external field are forced to 
obey Bfj-jt = BJxt = and Bf^^ = _B^j,^ = at the surface. It is thus immediately clear that we cannot match this kind of 
internal field with an external dipole field. One also has a condition on B*, as the exterior solution for a field which is regular 
at infinity is of the form 

B... = Y. f [-(^ + ^)^2yn + ei^doYn + 4>i^ sinerri (50) 

It is clear that if the f component must vanish at the surface, so must the 9 component. We thus have that all components of 
the mixed magnetic field must vanish at the surface. As we shall see in the following, this is true also in the case of a polytropic 
equation of state. In conclusion one can, therefore, match an interior mixed poloidal/toroidal dipole field only to a vanishing 
external field. It is also clear that the model does not allow for a purely toroidal field, since (I29|l is only compatible with the 
surface condition if c = 0. It is interesting to note how, as a consequence of the boundary conditions, we are restricted to a 
very limited class of magnetic fields: notably poloidal or mixed poloidal/toroidal fields that vanish outside the star. One may 
wonder if this restriction derives from having taken a somewhat pathological and simplistic model, the constant density one. 
To establish that this is not the case, let us move on to a more realistic case. 



4 POLYTROPIC MODELS 

Having understood the restrictions of the uniform density case, let us now allow for a non-uniform density distribution. It is 
natural to consider an n = 1 polytrope, as this model has many features of a realistic neutron star model. It also permits an 
analytic treatment. The allowed magnetic field must still satisfy equation (|4)), but now the role of p is no longer trivial as we 
are considering a density profile 

Pir) = Pc— sin (^) (51) 
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4.1 Poloidal field 

Following Monaghan 1 1996t ). we solve equation (|10p . imposing regularity at the centre and matching to an external dipoli 
The solution for the stream function is then 

which leads to a field of form 



(52) 



Br 

Be 



Bs cos 6 I ^3 



.3 _ ^2R^ _ .V^) sin (^) + 6 R\. cos (^)] ^ 



Bs sin 6 
'~2 



(53) 
(54) 



where Bs represents the conventional dipole field strength at the surface. 



4.2 Mixed toroidal and poloidal fields 

Let us now assume a mixed poloidal and toroidal configuration. We w ill once again loo k for dipolar solutions and take a 
stream function of the form S(r, 0) = A{r) so? 9. Furthermore, following Roxburgbl 1 19661) we shall define (3 to be 



(55) 



This is obviously a particular choice. However, it is natural since it leads to a separable equation which means that the problem 
can be treated more or less analytically. By varying the parameter A, which describes the relative strength of the toroidal 
part of the field, one can hope to get some insight into the how the two field components interact in the mixed problem. 
The magnetic field thus takes the form 

'2Acos6 -A' sin 6* AAsin^ 



B 



Then equation pO|) can be written 



A 



d_ 

dr 



1 f2A (fA 

pj.2 \^ J.2 ^j.2 



rR 



\'A 







(56) 



(57) 



pR'^r^ 

As before, we want to find a solution such that the field remains finite at the centre, cf. (|2ip . and that all the components of 
the magnetic field be continuous at the surface. As the toroidal field must vanish at the surface this forces the condition 



It must also be the case that 



A = at r = R 



A' = 



(58) 



(59) 

which derives from the condition that all components of the field must vanish at the surface. We can thus only consider fields 
that vanish outside the star. 

Solving the above equations for the density profile of an n = 1 polytrope, we find that the stream function takes the form 



Air) 



''"'^\ 2^3 ^2^2 /'"''"\2 >2 in f rX 

sm I —J r X cos A — sm J ^ ^ sm A — A sm ^ cos Att — 2 A Rtv cos I I 



-2 A Rtv r cos ( — 
R 



-2 sin 



rTr\ 

'rJ 



2 A sin ( — ^ R^ cos A + sin ^— ) sin Att^ + 2 Rtv r cos 
rX 



sin A 



R sin X + 2R TV sin ( — 
R 



rF{X) 



where Bk parametrises the strength of the magnetic field. We have defined 

F{X) = [A cos A -sin A] (A^ - Tv'^f 
Imposing the boundary conditions we find an eigenvalue relation for A. Permissible fields must be such that 

^ A(^2_^2) 

tan A = - ^ ' 



(60) 
(61) 
(62) 



7r2 - 3A2 

By taking higher values of A we can increase the relative strength of the toroidal part of the field compared to the poloidal 
part. Solving the trancendental equation numerically we find that the first three eigenvalues are 



Ai = 7.420 
Aa = 10.706 
A3 = 13.917 



(63) 
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These results agree with the values given in Roxburgh ( 19661 ) with an accuracy of « 0.1%. The key point is that the toroidal 
part of the field is not freely specifiable in this model. If we prescribe the strength of the poloidal field component, then we 
can find solutions with increasing toroidal fields, but these are discrete. 



4.3 Purely toroidal fields 

In the case of a purely toroidal field equation Q takes the form 

W- ( ^iB.rsine) - ^ ( ^(B.rsine) (64) 
or \prsmO J oO oO \prsmO J or 

The boundary conditions we have to impose are that the field vanish both at the centre and at the surface of the star, as 
before. We can take a solution of the simple form 

B ( TTT \ 

= — sin — sin (65) 

TV \ R / 

which will satisfy the boundary conditions since p = at the surface, as is the case for the polytropic model we are considering. 
Note that, in contrast to the uniform density problem, the polytropic equation of state allows for a purely toroidal field. 



4.4 Field confined to the neutron star "crust" 



It is often assumed that, if the core of a neutron star is a type I superconductor t he magnetic field will be expelled from 
the core and confined to the crust, see for example ^Bonazz ola fc Gourgoullion |l99^). To discuss this situation one should in 



principle consider the full equations of hydrostatic equilibrium, including the elastic terms. However, in order to investigate 
this problem, we will consider the case of a fluid with a magnetic field confined to a region close to the surface. This model 
would be relevant provided that the deformed shape represents the relaxed configuration of the crust. Then there is no strain, 
and thus no elastic terms in the equilibrium equations. We shall consider the same mixed field as in the previous section, 
i.e. of the form (|56|) . We shall, however consider the field to vanish inside a certain radius rt (which can be considered to 
be the base of the crust). The boundary conditions for the third order differential equation (|57|) at the surface thus remain 
A{R) = A'{R) — together with a condition which comes from imposing continuity of the Br component of the magnetic 
field, i.e. imposing 

A{rt) = (66) 

We again get an eigenvalue problem for the parameter A, allowing us to calculate the permitted ratios of toroidal to poloidal 
field strength. It should be noticed that we are not imposing continuity of the tractions at the inner boundary. In fact, there 
is a discontinuity in the B^ component of the field at rj, which will lead to currents at the crust/core interface. In order to 
simplify the calculation we also take the core to be unperturbed, and simply impose continuity of the perturbation in the 
gravitational potential 5$ and of its derivative 5^' . 

Before moving on, it is worth remarking that this simple model can only be seen as a rough representation of the problem 
with a superconducting core. There are a number of issues that one ought to worry about, and which we are not considering 
here. For example, the complete flux expulsion from a Type I superconducting core may be a severe oversimplification. Is it 
really to be expected that the magnetic flux is expelled from the entire core rather than (say) bunched up into macroscopic 
regions? Moreover, if the core forms a type II superconductor (as is usually expected) then the magnetic flux is carried by 
fluxtubes. This problem is signiflcantly different from our idealized model, and requires a separate analysis. 



5 MAGNETIC DEFORMATIONS 

So far we have discussed how the requirement that the magnetic field to be consistent with the fluid conflguration greatly 
restricts the model parameter space. In particular, it does not allow us to freely specify the toroidal component of the field in 
the mixed case. Also, for the models that we have considered, one can only consider fields which vanish in the exterior. This 
is clearly not physical, as one would expect the exterior field to be prevalently dipolar far from the star. This problem could 
easily be fixed by allowing surface currents (which one may expect as there are strong electric fields at the stellar surface). 
Nevertheless, given that we do not have a physical model for such currents we will not introduce them here. 

Having determined some magnetic field configurations which are consistent with a chosen stellar interior, we can turn 



our at tention to solving equations ([T| in order to obtain the new equilibrium shape of the star. Inspired by the work of ISaio 



( 19821 ) and our own recent work on crustal deformations (|Haskell et al.ll2006l ). we shall define a new variable 



x{r,e)^r[l + e{r)Pi{e)] (67) 
where r is the standard radial variable and Pi is a Legendre polynomial representing the deformation. Note that we are 
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considering only one polynomial here, but the formalism could easily be extended to a sum of Legendre polynomials. The 
perturbed surface thus takes the form 

Xs{R,9) ^ R[l + e{R)Pi{9)] (68) 

where R is the radius of the unperturbed (spherical) star. We shall also assume that the pressure, density and gravitational 
potential take the form 

iP{r,6) = iP{r) + 5tp'{r)Pi (69) 

where ip is the background quantity and Sip' is a small perturbation of O(B^). From now on we shall, unless there is a risk of 
confusion, write Sip instead of SipK Equations ([TJ thus take the form 

+ (70) 
dr dr dr J 4tv 

iSp + pS,>)vp = l^^2L^l2L^ (71) 

and must be solved together with the perturbed Poisson equation 

d^S^ 2dS^ 6 , , 

+ --7^-IT'5* = 4^5p (72) 



dr^ r dr r 
Let us first of all tackle the case of an incompressible star 



5.1 Deformations of incompressible stars 

In the case of a uniform density star we consider the field in equation (|26p . This gives us for the Lorentz force 
2 f bV f R^r 2r^ 



as well as, obviously, = 0. In the case of an incompressible star, there can be no / = deformation, as this would not 
conserve the volume and therefore not conserve the mass of the star. We shall thus only consider the case of / = 2. For this 
case 5$ inside the star takes the form 

Sii> ^ --TTGpe(Ry (75) 
5 

Inserting this solution into the 6 component of (|71|) and evaluating at the surface gives 

5p{R) = j^7vGpe{R)R:' - ^ (76) 

If we now remember that 



we have at the surface 



By using the background pressure 



5p{r) — 5p{x) — £{r)r — {r) (77) 
dr 

Sp{R) ^0-e{R)R^iR) (78) 
dr 

p = 2nGp'^{R^ -r^), (79) 

which agrees with the result of GoosensI ( 1972 ) . 

In order to assess the relevance of this deformation for gravitational-wave emission we calculate the ellipticity, which is 
defined as 

Here Iq is the moment of inertia of the spherical star, while Ijk is the inertia tensor 



equation (I76[) gives us 



(81) 



Ijk = / p(r){r'^Sju - XjXk)dV. (82) 
Jv 
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For a constant density star this leads to 

3 _IB^R^^ R Y f M y f B Y 

' - "2^-T8GaF~^° [Wk^J [Tam^J [Ta^J ^^^> 

where B is the energy averaged magnetic field. The e llipticity is po sitive, so the star is oblate, as expected. As a sanity check 
of this result we can compare to the estimate used by Cutler ( 2002 ^1 . For a field of 10^^ G he assumes that e ~ 1.6 x 10~®. Our 



more detailed calculation has led to a result that is almost a factor of two smaller than this. Nevertheless, the two estimates 
are reasonably close. 

5.2 n=l polytrope with a poloidal field 

Let us now consider the case of a star with an n = 1 polytropic equation of state. The magnetic field appropriate for this case 
is that in equation (I54p . for which the Lorentz force has non-trivial I — 2 component^ 

Lr = -I [2r\' + - mr\-) sin (^) + (3r^.^ - m\.) cos (!g)] ^ r {^) ^^(^) ^ 



■ (7r2 - 6) . 

-i,' [rV - - (II) + 6K'„co, (5)] .in (^) ^ ,85, 



2 

s 

■ (7r2 ^) _ 

Having worked out the Lorentz force we can solve the perturbed hydrostatic equilibrium equations (|71|l and the Poisson 
equation (|72[). with the condition that 5^ and S^' be regular at the centre and match the exterior solution at the surface. 
From equation (|7ip we then obtain 

B^TT 



Sp^- pTr" + {SRr^' - 6R\'') sin (^) + 6 7V^R\ cos (^) + 2pJ^ (r) (^-^^ r(7r' - 6)'] 



(86) 



6)2Gpci?* 

LVcosf — ) -2rV-(3rVi?-36i?^ + 12 7rVii^)sinf — ) -36^ri?^cosf— )1 (87) 



Inserting this into the right hand side of equation (|72|l we can solve for 5^ 



5$ 

We can then evaluate the distortion at the surface 



Sp dp TT^ (tt^ - 24) 



2 



16i?2p^2(;j(7r2 -6)^ 

Finally we obtain the ellipticity 

■K^R (37r2 - 32) R'^Bs'^ 



(89) 



(7r2 -6)^ GpcM2 

Scaling to canonical neutron star values for the different parameters, and introducing the volume averaged field as in the 
constant density case (here B — 0.54Bs), we have 

— »- (is^^y (i^)"' (bSs)' 

In other words, for this model configuration the magnetic field induced ellipticity is almost two order of magnitude greater 
than in our constant density model. This is perhaps a useful indication of how much the result can vary for supposedly 
"similar" field configurations if we change the equation of state. Of course, in making this statement one must keep in mind 
that the two field configurations really are different. 



5.3 n=l polytrope with a toroidal field 

We next consider the case of a purely toroidal magnetic field in an star with an n = 1 polytropic equation of state. For the 
field given in (|65p the Lorentz force has I = 2 components, 

U = |^[^sin(!J)+.rcos(^)]sin(^)P2(e) 



2 In the general case there will be both a radial and a quadrupole deformation. Since we are primarily interested in the gravitational- 
wave aspects we ignore the I = deformation in the discussion. This contribution would be important if one wanted to (say) study the 
oscillations of deformed magnetic stars. The required calculation is provided in Appendix B. 
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From equation (|71|l we can then write 

^P = ~ 247?3pcG7r 
which, inserted on the right hand side of Poisson's equation (|72p gives 



5* = (I) ' (rV^ - 15rii*) cos (^) + 5R' (R 



36 V-R 

We can then calculate the distortion at the surface 



It follows that the ellipticity is given by 



r TT sin 



7rr \] 1 



1 TT^ 



15 



144 n^R^Gpc 



1 (tt^ - 15) 7?" 



97r (7r2 - 6) Gp^ 

or using the average field, which in this case corresponds to B ~ — 0.17-B, 



-5 X 10" 



R 



10 km 



M 



1.4 M© 



B 



1012G 



(92) 
(93) 
(94) 
(95) 

(96) 



The deformation is now prolate, and notably about one order of magnit u de sm aller than in the case of a poloidal field (for 
the same B). We can compare this result to the estimate used by Cutler ( 20021 ) . For a 10^^ G toroidal field he assumes that 
e « —1.6 X 10~^ (the same magnitude as in the poloidal case). The deformation we have determined is about a factor of three 
larger. 



5.4 General deformations 



Having considered some particular examples, let us now present the formalism for the case of a more general field configuration 
and equation of state. Assume that the magnetic field takes the form of equation (|56|l . for which the Lorentz force is 



Le 



dA sin I 
dr 

2A 



(A^ 



2)A + 



d'A 
dr^ 



■ cos U sm t 



,,2 2 N d A 



(97) 



Focusing our attention on the I — 2 components, we have (see Appendix B for a discussion of the associated radial deformation) 



Lr{l = 2) = 



dA 2P2{e) 
'dr 3r4 
2A dP- 



\'Ar' + 



Lb{1 = 2) = =—-^[X^Ar^ + — ^2A 



3r5 de 



d^A 
dr^ 
d^A 



2A 



(98) 



From the angular part of the perturbed hydrostatic equilibrium equations ^ we obtain (for the 1 = 2 case we are considering) 



5p = — p5<E> + rLg 
2 A 



d'A 



dr"^ 



which substituted back into the radial part of the equations leads to 



5p 



1 



3r5 4£ 

dr 



dp 



5<i,r'' + r'^A (2r^A^ - 4) + AA' (4 - r^A^) + 2r^ A 



ArA- 



.(fA 



dr dr dr^ dr'^ 

This now allows us to compute the source term of the perturbed Poisson equation (|72p . which, for / = 2, reads 



d^5^ 2 dS$ 



dr^ ^ r dr 



<5$ 



6 4ttG (dp 

-7+ ' 



54* \dr 



AtvG 

dr 



r-^yl(2r^A' - 4) + 4yl'(4 - r^A^) + 2r^A 
dr 



(£A 

dr^ 



'4r^A 



d'A 
dr^ 



(99) 



(100) 



(101) 



The boundary conditions for this equation remain regularity at the centre and continuity of 5$ and its derivative at the 
surface. For I = 2 the exterior gravitational potential has the form 

1 



(102) 



Having computed 5$ in the interior we can obtain Sp from equation (|100|l and thus compute the ellipticity e and the 
deformation at the surface 

Sp{R) 



eiR) 



(103) 
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A 


e 




Bt 


7.420 


-6.3 X 10" 


13 


3.0 X 10^2 


10.706 


-2.8 X 10" 


12 


4.2 X 10^2 


23.433 


-4.4 X 10" 


11 


8.1 X 10^2 


80.073 


-2.3 X 10" 


10 


2.7 X IQis 


105.215 


-4.4 X 10" 


10 


3.5 X 10^3 


183.994 


-3.7 X 10" 


8 


9.3 X IQis 


1000.59 


-1.3 X 10- 


6 


4.8 X IQi"* 



Table 1. The ellipticity e for a sample of eigenvalues A. The results correspond to the parameters M = 1.4M0, i? = 10 km and we take 
the energy averaged amplitude of the poloidal part of the field to be 10^^ G. The energy averaged toroidal field is given in each case. As 
can be seen, the star starts off with a small deformation as the effects of the poloidal and toroidal components cancel each other, and 
becomes more and more prolate as the toroidal component of the field grows. Note that we have a vanishing exterior field in all cases. 



A e Bp Bt 

103.60 -1.0 X 10"^ 6.2 X lOi"* 1.2 x lO^s 
496.39 -5.5 x 10"^ 5.8 x IQ^^ 1.3 x lO^^ 
1008.46 -1.7 X 10"5 5.6 x 1.4 x lO^^ 

Table 2. The ellipticity e for various eigenvalues A. The results correspond to the parameters M = 1.4Mq, i? = 10 km. We take the field 
to have the same energy as that of a field extended to the whole star with A = 1000.59. The energy averaged values of the poloidal and 
toroidal magnetic field are also given. We should thus compare all results for the ellipticity in this table with the one given by A = 1000.59 
in Table [T] The field is confined to a region close to the surface beginning at a radius r;, = 9 X 10^ cm, thus roughly corresponding to a 
1 km thick crust. The core is taken to be spherical. As can be seen for the listed results, the deformations are larger than in the case of 
the field extending to the whole star, by up to an order of magnitude. 



As an example we carry out the calculation for a polytrope and a mixed poloidal/toroidal field. Then we can use the 
stream function from equation (|60p to determine the magnetic field and the Lorentz force. The result will obviously depend 
on which of the eigenvalues we choose for A. This way we can vary the relative strength of the toroidal part of the field. A 
sample of the obtained results is listed in Table [T] As expected, for low values of A the effects of the poloidal and toroidal 
field work against each other and the deformation is comparatively small. It is possibly surprising that the cancellation is so 
efficient. If one were to, for example, add the deformations deduced for a pure poloidal field (|90[) to that for a purely toroidal 
field (|96|) using the appropriate values for B then one would predict that the deformation should be e ~ 10~^^ for the first 
eigenvalue in the Table. Furthermore, one would have expected the deformation to be prolate. We see that for the mixed field 
we are using the true deformation is almost two orders of magnitude smaller. It is also toroidal. This illustrates that it will be 
important to use both realistic field configurations and equations of state in estimated magnetic neutron star deformations. 
The results in Table [1] show that the star is always prolate, and as A grows the toroidal field becomes dominant and the star 
becomes more and more prolate. Once a toroidal component is introduced, the poloidal component of the field is, thus, never 
strong enough to make the star oblate for this model. 



5.5 n=l polytrope, field confined to the crust 

Finally, let us consider the case when the field is confined to a region close to the surface. As we have already discussed, 
we shall call this "crust", even though we are neglecting its elastic properties. We consider the core to be unperturbed, and 
impose the condition (|66|l at the inner boundary. Then the deformation is confined to the crust and, as the density is low 
in this region, the quadrupole moment, and thus the ellipticity, will in general be much smaller than if we were considering 
the deformations of the whole star. However, it is possible that the field is confined to the crust after being expelled from a 
superconducting core. In this case a large number of field lines are squeezed into a small region close to the surface of the star 
and would produce a strong field and large deformations. This result can be seen in table [2 where we have calculated the 
deformation of the star by assuming that the total magnetic energy in the crust is equal to that of the field extended to the 
whole star, as calculated in the previous section. As can be seen the ellipticity is larger than in t he case of a field extending 



throughout the star, by up to more than an order of magnitude. This agrees with what was found in lBonazzola fc Gourgoulhon 



(|l996), where the authors analysed the deformations due to a poloidal field confined to a thin shell close to the surface of a 
neutron star. This motivates a more detailed investigation into the problem of superconducting neutron star cores. 



6 CONCLUDING DISCUSSION 

We have presented a scheme for calculating the magnetic deformations of a neutron star. We have considered the particular 
case of a dipolar field in a non-rotating star, but the formalism is readily extended to other cases. The formalism can, for 
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example, be adapted to the case of slow rotation if we take the magnetic axis to be aligned with the spin axis of the star 
(see Appendix A) . We have considered the case of a purely poloidal field, the effect of which is to make the star oblate, that 
of a purely toroidal field which makes the star prolate, and the case of a mixed poloidal and toroidal field, which makes the 
star more and more prolate as the strength of the toroidal component is increased. We have seen that the condition that the 
magnetic field be consistent with the background equation of state is very restrictive. For the present set of model problems, 
we are forced to consider only fields that vanish outside the star. This is clearly not physical as one would expect the field to 
be prevalently dipolar far from the star. This problem could be solved by introducing surface currents. Since we do not have 
a physical model for such currents we have not considered this possibility in detail. Nevertheless, our formalism is sufficiently 
general that it could be extended to any field configuration. 

It seems appropriate to conclude this investigation by discussing the relevance of our results for gravitational-wave 
observations. A rotating magnetic neutron star becomes interesting from the gravitational-wave point of view when the 
magnetic axis is not aligned with the spin axis. Then the deformation is no longer be axisymmetric and we have a time 
varying quadrupole. This is in contrast with the typically much larger (see Appendix A) rotational deformation which is 
always axisymmetric and cannot radiate gravitational waves. The framework we have described can, in principle, be used 
with any equation of state or field configuration for the neutron star. It allows us to calculate the ellipticity, and thus 
the quadrupole, in order to make quantitative gravitational-wave estimates. This is a useful exercise since we can compare 
our results to observational upper limits on the neutron star ellipticity, and perhaps constrain the parameters of proposed 
theoretical models. T o see how this work s out in the present case, let us consider the recent upper limits on pulsar signals 
from the LIGO effort ( Abbott et al.l[2007 ). The strongest current constraint comes from PSR J2124-3358, a millisecond pulsar 



spinning at 202.8 Hz. For this system, LIGO obtains the constraint e < 8.5 x 10~^. This is impressive, because it restricts the 
maximal neutron star mountain to a fraction of a centimetre. The question is, how much does this constrain the magnetic 
field? From (|90p and (|96[) we easily find that we must have 

J 10^* G for a purely poloidal field (104) 
j 4 X 10" G for a purely toroidal field ^ 

(taking radius and mass to have the canonical values) . Is this constraint telling us anything interesting? The answer is probably 
no. The pulsars currently considered by LIGO are almost exclusively fast spinning, taken from the millisecond spin period 
sample. For these neutron stars the exterior dipole field inferred from the spindown is weak. In the case of PSR J2124-3358 
the standard estimate would suggest that the exterior field is about 3 x 10* G. To understand the relevance of this we must 
first recall that our model stars have a vanishing exterior field. Thus there is something missing in the model. However, it is 
thought that the interior field could be significantly stronger than that in the exterior. It could well be that our simple models 
provide useful representations of this interior field, and what is missing is a relatively small part that is irrelevant when it 
comes to deforming the star. Nevertheless, it does not seem reasonable to suggest that the interior field can be as much as six 
orders of magnitude stronger than the exterior one (for reasons of stability etcetera). Hence, our results indicate that LIGO 
has not yet reached the sensitivity where it constrains any reasonable theory. Given this, one may ask to what extent future 
observations will improve on the current results. It is relatively straightforward to estimate how well one can hope do to if we 
recall that (using matched filtering) the effective gravitational-wave amplitude increases as the square-root of the number of 
observed wave cycles. For a spinning neutron star, where the spin frequency is essentially fixe d during the observa tion, this 



1/2 A ^ 

means that the signal-to- noise ratio scales with the observation time as t^[^^. The data used bv lAbbott et al.l (|2007l ) is based 



on observations lasting the order of three weeks. If one could extend this to one full year, one would gain a factor of 4 or so 
in amplitude. This translates into an improvement of a factor of 2 in the constraint on the magnetic field, not enough to test 
a realistic theory model. In order to do much better one would have to improve the detector sensitivity. With the advanced 
LIGO upgrade, the sensitivity will increase by about one order of magnitude. With a one year observation we then get a 
further vTO improvement in the constraint on the magnetic field. Advanced detectors will also have the ability to narrow-band 
and focus on a particular frequency. Suppose one were to do that and get (say) an additional order of magnitude improvement 
of sensitivity at the particular frequency for a given pulsar. This would still only improve the constraints in (|104p by a factor 
of about 2 X vTO x \/T0 — 20. These rough estimates suggest that we should probably not expect to detect gravitational 
waves from magnetically deformed millisecond pulsars. There are obvious caveats to this, in particular, associated with the 
superconductivity of the core. If all the magnetic flux is confined to the crust then the deformation may be larger, cf. the 
results in Tabled This problem is clearly worth further investigation. 

The situation is not quite as pessimistic if o ne considers slower spinning pulsars. Consider for example the current LIGO 



upper limit of e < 2 x 10 for the Crab pulsar (jAbbott et al.ll2007l \ For the models we have considered we must then have 



, 5 X 10^^ G for a purely poloidal field 

' 2 X 10^^ G for a purely toroidal field ^ ^ 

At first sight, these limits may seem less interesting than those from the faster spinning system. However, we need to remember 
that young pulsars, like the Crab, are thought to have much stronger magnetic fields. From the spin-down one would usually 
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estimate that the exterior dipole field is 4 x 10^^ G for the Crab. Hence, the current observed upper limit is only a factor of 
about a factor of 1000 away from the expected field. One should be able to improve this by at least an order of magnitude 
with future detectors (the current LIGO detectors may be improved at lower frequencies, Advanced LIGO will do better and 
VIRGO is expected to have good performance in the low-frequency regime) and a full year of observation. This suggests that 
gravitational-wave observations may in the future be able to test theoretical models where the interior magnetic field is more 
than two orders of magnitude stronger than the exterior field. 

The challenge now is to improve on the rough estimates we have presente d in this paper. Futur e work needs to consider 
the deformation for realistic field configurations such as those worked out by iBraithwaite fc Spruit We also need to 

understand the role of superconductivity better. There are a number of interesting questions associated with accreting systems. 
These involve the bur ial of the field by accre ted matter, leading to a relatively weak exterior field, and potential accretion 
induced asymmetries I Melatos fc Paviig 2005h . Finally, it may also be interesting to consider magnetars. After all, they are 
expected to have fields as strong as the contraints given in (|104p . Of course, they are also spinning very slowly which places 
them outside the bandwidth of any groundbased gravitational-wave detector. It is still possible that a newly born magnetar 
may occasionally pass through the detection window. Such an event should be observable from within the galaxy, and perhaps 
beyond. 
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APPENDIX A: GENERAL DEFORMATIONS: MAGNETIC FIELDS AND ROTATION 



The framework discussed in the main body of the paper can easily be extended to the case where the deformation is due not 
only to a magnetic field, but also to the star's rotation. We can write then the equations of motion as 

L 

47rp 

where 



p Aixp 



(Al) 



= $ - l/2nV%in^6' 



(A2) 

$ is the gravitational potential and the constant rotation rate. As the new term due to rotation is written as the gradient 
of a scalar function, it's curl will still vanish and the compatibility condition for the magnetic field in equation @ remains 
the same. We can thus still use the magnetic field from equation (|56p . The equations of hydrostatic equilibrium, for the I = 2 
case, now take the form 

dSp dS$ , 2 . 

— + p—— + dp— + -pQ. r^Lr 
dr dr dr 6 



Sp + pS^ + 



pQf-r 



rLe 



(A3) 



where Lr and Lg are given in equation (|98|) . We can proceed as in the previous analysis and use equations HA3|) to obtain Sp 
as a function of 5$, and ft. This allows us to write the perturbed Poisson equation (for 1 = 2) as: 



d'^5^ 2 dS$ 
dr'^ r dr 



5$ 



_6_ 

^2 ^T:,d± 



4-kG fdp 
dr 



r—A(2r 
dr 



4) +4yl^(4-r^A^) 



2r-'A 



fA 
dr^ 



^^2^d A ^ dp 
dr^ dr 



il r 



(A4) 



and thus obtain the surface deformation e and the ellipticity e as defined in equations (|102|) and (|103[) . 

As an example let us take an n = 1 polytrope with the purely poloidal field from equation (|54[) . The equations of 
hydrostatic equilibrium now include the rotation term as in equation (|A3|I . where the Lorentz force is that of equation (|85p . 
We can thus apply exactly the same procedure as previously and solve for the perturbed quantities 5$, 5p and 5p. This allows 
us to calculate the deformation at the surface 

_ 7r^B^(24-7r^) 5 »^ 
^ ^ " 167?2p?G(7r2 - 6)2 " 4^G^ 



(A5) 
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Working out the corresponding elhpticity we find 



R V f M X"* / B V , . .^-T ^ ^ 



^ = [Tol^) [Tam^) [m^) (^^^ 

where Qir = |\/7rGp is the breakup frequency (p is the average density of the non-rotating model). For our values this 
frequency would be / « 1250 Hz, which corresponds to a period P « 0.8 ms. We see that, as expected, the rotational 
term completely dominates the magnetic term. Rotational deformations will, however, always be axisymmetric. The magnetic 
deformations can lead to a quadrupolar deformation and thus to gravitational wave emission only if the magnetic axis is 
inclined with respect to the spin axis. 



APPENDIX B: RADIAL DEFORMATIONS 



As discussed in Section 5, it may be relevant to work out also the radial deformation due to the presence of the magnetic field. 
For general fluid configurations, any quadrupole deformation is likely to be accompanied by an Z = component. In order to 
determine the radial deformation we can assume that the new radial variable x defined in equation (|67|l labels the deformed 
gravitational equipotential surfaces. This mean that we impose that 

Va; X V$ = (Bl) 

which, if we write 

e{x, 9) = Do{x) + D2{x)P2{e) (B2) 

leads to the condition 



D2 = 



(B3) 



This expression allows us to calculate the quant i ty D2 throughout the star and can be of use, for example, when discussing 
oscillations on a deformed background, cf. I Said (119821 1. Writing out the equations that need to be solved in the I — case, 
the perturbed Euler equation and Poisson equation, we have 



dr 

-7-(5i/'o 
dr 

■dr^P' 



2 

AnGSpo tpo 



£_ 

P 



Lo{B) - pifjo - 



dr\p) ^ ^ dr 



(B4) 



where ipo is introduced in order to obtain a first order system. We have also assumed a linearised barotropic equation of state 

(B5) 



5P Sp 



and Lq{B) is the Z = part of the Lorentz force, which can be calculated once the magnetic field has been specified. These 
equations can now be integrated, imposing regularity at the centre of the star and imposing that at the surface Spo = and 
matching (5$o and S^'q to an exterior solution of the form J^™' oc 1/r. This will allow us to calculate 

-(5$o 



Do 



4n- 



(B6) 
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